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ABSTRACT 



The evolution of the frequency-directional wave spectrum, E(f,d), across the inner 
continental shelf and beach was examined with measurements collected at the U.S. Army Corps 
of Engineer’s Field Research Facility during the recent SandyDuck experiment. Arrays of bottom 
pressure sensors were deployed on the shelf in 20 m depth and on the beach in depths ranging 
from 2 - 5 m. These arrays were complemented by a directional wave buoy in 20 m depth and 
an array of pressure sensors in 8 m depth maintained by the U.S. Army Corps of Engineers. A 
preliminary analysis of these data is presented here focused on four case studies that illustrate the 
observed wave shoaling evolution in both non-breaking and breaking conditions. Estimates of 
E(f,6) extracted from array cross-spectra at six cross-shore locations are compared to predictions 
of linear refraction theory. The present observations support conclusions from previous studies 
that the cross-shore evolution of dominant wave propagation direction is well described by linear 
refraction theory. Observations of harmonic peak development at directions aligned with the 
dominant waves are consistent with theoretical wave-wave interaction rules and previous 
observations. In both non-breaking and breaking conditions, the observed E(f,0) are directionally 
broader than predicted. In contrast to previous observations on a barred beach, the present 
observations on a planar beach do not show a dramatic broadening of directional wave spectra 
in the surf zone. 
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I. INTRODUCTION 



Accurate predictions of the nearshore transformation of ocean waves are vital to 
naval littoral operations (e.g., mine warfare, amphibious landings, swimmer operations) 
and civilian coastal erosion studies. As waves propagate from the open ocean onto 
shallow coastal shelves and ultimately break on beaches, linear and nonlinear effects 
change their amplitudes and directional characteristics. The present study is focused on 
the still poorly understood transformation of the directional properties of waves 
propagating across the inner shelf and beach. 

Linear shoaling of waves is well understood. Reduction of the group speed in 
shallow water generally causes an amplification of wave amplitudes on beaches. Waves 
arriving at an oblique angle relative to the shore are refracted toward normal incidence 
while for large incidence angles their amplitudes are reduced (e.g.. Kinsman, 1965). 

Close to shore, strong nonlinearities and wave breaking affect the frequency 
spectra and directional properties of waves. Nonlinear interaction between two primary 
wave components with frequencies /j and ^ and vector wavenumber ki and kz excite a 
secondary wave component with the sum frequency (/1+/2) wavenumber (ki+kz). If 
the incident wave spectrum is narrow (i.e.,/j~/2» ki^kz), these secondary waves are 
harmonic components with approximately twice the frequency of the incident waves and 
propagate in approximately the same direction as the incident waves (Freilich et al. 1990; 
Berbers et al. 1998). Subsequent interactions between the primary waves and the newly 
formed secondary waves excite higher harmonic components (3/j, 4/j, ...) all 
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approximately directionally aligned with the incident waves. 

Field observations of the shoaling evolution of directional wave spectra are scarce. 
Freilich et al. (1990) presented high resolution directional spectrum estimates obtained 
from two arrays of pressure sensors located weU outside the surf zone on a nearly planar, 
natural beach. These observations show the expected refraction toward normal incidence 
of incident swells and energy transfers to directionally aligned harmonics. In bimodal 
spectra, with two distinct wave systems arriving at the beach from different incidence 
angles, the theoretically expected excitation of waves at the sum frequency and vector 
wavenumber was observed. Berbers et al. (1998) presented estimates of mean directions 
and directional spreading parameters extracted from cross-shore transects of co-located 
pressure sensors and current meters on barred and non-barred beaches. These 
observations show that although the mean wave direction does not appear to be strongly 
influenced by wave breaking, the directional spreading of waves increases significantly 
in the surf zone on both barred and non-barred beaches. 

In the present study, a preliminary analysis is presented of a more extensive 
shoaling wave data set acquired during the SandyDuck experiment near Duck, NC. Six 
high resolution arrays of bottom pressme sensors were deployed for fom months along 
a transect extending from the inner shelf through the surf zone. These measurements 
provide unique detailed estimates of the shoaling and breaking transformation of 
directional wave spectra. The field experiment and data analysis techniques are described 
in chapter II. Four case studies spanning a range of conditions (i.e., nonbreaking and 
breaking waves, narrow swell spectra, and broad wind seas) are described in chapter III. 
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The results are summarized in chapter IV. 
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n. FIELD EXPERIMENT AND DATA ANALYSIS 



Detailed measurements of the evolution of frequency-directional wave spectra 
across the inner continental shelf and beach were collected during the SandyDuck 
experiment at the U.S. Army Corps of Engineer’s Field Research Facility located near 
Duck, NC, on a straight barrier island. This site is fully exposed to the Atlantic Ocean 
and features an 80 km wide, shallow (20 - 50 m depths) shelf. Arrays of pressure sensors 
were deployed on the inner shelf in 20 m depth (5 km from shore) and on the beach in 
depths ranging from 2 - 5 m (Figure 1). The experiment was conducted in collaboration 
with Dr. Steve Elgar (Washington State University) and Dr. R. T. Guza (Scripps 
Institution of Oceanography). High-quahty wave data were collected at all instrument 
locations during a fom month period (August - November 1997). 

The 20 m depth array consisted of nine pressure sensors arranged in a 500 m 
equilateral triangle (Figure 2, bottom). The pressure sensors are battery powered, internal 
recording instruments mounted rigidly near the sea floor inside anchors. The autonomous 
instruments were synchronized to within 10 ms using accurate temperature compensated 
oscillators and thus can be used as a coherent array for sea and swell waves. The sample 
frequency was 2 Hz. Instruments were accurately positioned using differential GPS, an 
underwater acoustic navigation system, and direct underwater tape measurements. Sensor 
spacing uncertainties range from about 0.4 m for the shortest array legs to 2.0 m for the 
longest array legs. The array geometry allowed for omni-directional wave measurements. 
Sensor spacings ranging from 20 m to 300 m were chosen to accurately resolve swell and 
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sea waves in the frequency range 0.04 Hz - 0.2 Hz. Waves with frequencies higher than 
0.2 Hz are strongly attenuated at the sea floor, and thus not accurately measured by the 
pressure sensors (e.g., Herbers and Guza, 1991). A Datawell Directional Waverider buoy 
was deployed within the array perimeter to provide measurements of the high frequency 
part (0.2 Hz - 0.5 Hz) of the wave spectrum. This surface-following buoy measures 3- 
component accelerations that yield wave frequency spectra and low-resolution directional 
information (e.g.. Long, 1980). Buoy data sampled with a frequency of 1.28 Hz was 
transmitted continuously via an HF radio link to a receiver on shore. 

The nearshore instrumentation included four arrays of 10 or 1 1 bottom-mounted 
pressure sensors, each with an alongshore aperture of 210 m and a cross-shore apertme 
of 35 m (Figure 2, top). The arrays were positioned at offshore distances (relative to 
F.R.F. coordinates) of 210, 260, 375, and 500 m in nominal depths of 2.9, 3.7, 3.9, and 
5.1 m, respectively. The nearshore arrays also included sonar altimeters that provided 
accurate, continuous water depth measurements (Gallagher et al. 1996). All instruments 
were cabled to a central data acquisition center onshore. The sample rate of the pressure 
sensors was 2 Hz. 

Additional directional wave data were available from a permanent array of 
pressure sensors (maintained by the U.S. Army Corps of Engineers Field Research 
Facility) located about 900 m from shore in 8 m depth (Figure 2, middle). The 15- 
element array with an aperture of 250 m (alongshore) x 120 m (cross-shore) and a 2 Hz 
sample rate resolves waves with frequencies up to approximately 0.3 Hz. 

Bathymetry at the site is characterized by a gently sloping (» 1:250), nearly planar 
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inner shelf (Figure 1) and a slightly steeper (-1:100) beach (Figure 3). Seaward of the 
20 m array, the shelf is characterized by a very gentle slope (-1:2500) with large (scales 
0(5 km)) ridge-like features that have vertical scales of 0(5 m). Refraction computations 
show that the effects of the benign shelf topography on wave transformation across the 
shelf are significant, but do not produce strong local gradients in wave energy 
(Hendrickson, 1996). Daily bathymetric surveys of the nearshore region (e.g.. Figure 3) 
indicate that the beach changes over the course of the experiment were small. The sand 
bar located 300 m offshore (Figure 3) remained relatively stationary with its crest 
submerged about 5 m below mean sea level, and did not strongly affect the wave shoaling 
process. Close to shore, the morphology was more dynamic, with the development of 
transient bars and alongshore depth variations. These changes occurred well shoreward 
of the nearshore arrays, and thus did not influence this study. 

In the following analysis of wave shoaling evolution from 20 to 2.5 m depth, 
alongshore depth variations are neglected. The preliminary analysis presented here is 
focused on four 3-hour-long case studies that illustrate the observed wave shoaling 
evolution in both non-breaking and breaking conditions. Prior to the analysis, the data 
were tested for stationarity and spatial homogeneity using a variety of methods. Temporal 
changes in 2-minute average wave variances were used to ensure stationarity of wave 
conditions throughout the selected 3-hour-long record. • Alongshore homogeneity over the 
alongshore extent of the nearshore arrays was verified by intercomparing the measured 
spectra of each sensor in the array, as well as the coherence and phase spectra of 
redundant alongshore sensor pairs. For example, the southernmost sensors of the 
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nearshore arrays consistently recorded spectral density values an order of magnitude lower 
than other array sensors for waves arriving from large, oblique southerly directions. 
These discrepancies are likely caused by the close proximity of the U.S. Army Corps of 
Engineer’s research pier (e.g., wave blocking by pier pilings), and therefore data from the 
southernmost sensors was not used. 

Array cross-spectra were calculated based on Fourier transforms of overlapping 
1024 s segments. Merging 1 1 frequency bands resulted in estimates with a frequency 
resolution of 0.0107 Hz and 220 degrees of freedom. For each frequency, /, the 
corresponding directional distribution of wave energy, S(9‘/), was estimated from the 
array cross-spectra using the variational technique described in Herbers and Guza (1990). 
The direction 9 is defined relative to the local shoreline orientation with 9=0’ 
corresponding to normal incidence to the beach (waves arriving from 70" true N), and 9 
is positive (negative) for waves approaching the beach from northerly (southerly) 
directions. The method yields the smoothest non-negative distribution fimction that is 
statistically consistent with the observations. This roughness namirriization technique 
effectively suppresses spurious peaks in S(Q'f) estimates that might otherwise be 
misinterpreted as multimodal features of the directional wave spectrum. The estimates 
of S(9i^ were combined with an estimate of the surface elevation frequency spectrum E(f) 
(obtained by averaging the measured auto-spectra and applying a linear theory depth 
correction) to form the wave frequency-directional spectrum E(^,9) (= E(/)S(9/)). 

Estimates of E(/',9) were initially computed at all array locations using the full 
array geometry (i.e., both the alongshore and cross-shore legs). Spurious peaks and poor 
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convergence of 8(0^/) estimates at high frequencies observed at the nearshore and 8 m 
depth locations, suggest that nonlinear and/or wave breaking effects contributed significant 
errors. The estimation technique assumes a linear dispersion relation, and thus is sensitive 
to nonlinear changes in the wavenumber (e.g., amplitude dispersion and bound waves). 
Wave-breaking induced energy losses may reduce the coherence in the cross-shore array 
legs causing additional inversion problems. Frequency-directional spectra at the nearshore 
and 8 m depth locations were subsequently re-computed using only the alongshore legs 
of the arrays. These linear array estimates are insensitive to dispersion errors and cross- 
shore inhomogeneities, but cannot distinguish shoreward propagating waves (-9O°<0<9O°) 
from seaward propagating waves (9O°<0<27O°). Since the full array estimates already 
confirmed that swell reflections from shore were weak, the S(0i/) estimates obtained from 
the linear arrays were constrained to vanish at seaward propagation directions 
(9O°<0<27O°) (Herbers and Guza, 1990). The resulting estimates were nearly identical 
to the full array estimates in the energetic part of the spectrum and smooth at higher 
frequencies. In two cases, August 17 and October 19, estimates of E(^0) extracted from 
the 20 m depth array measurements also indicated spurious structure and convergence 
problems at high frequencies. The estimates were recomputed using a five element, linear 
sub-array. In the August 17 case, with waves arriving from southerly angles, the northern 
leg of the triangle array was used with the constraint S(0;^ = 0 at angles (-18O°<0<-15O°, 
3O°<0<18O°). Similarly, for the October 19 case with waves arriving from northerly 
directions, the southern leg of the triangle array was used with the constraint S(Q J) = 0 
at angles (-18O°<0<-3O°, 15O°<0<18O°). 
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Estimates of E(40) in 20 m depth over a wider frequency range (0.04 - 0.5 Hz) 
were obtained from the directional wave buoy using the Maximum Entropy Method 
described by Lygre and Krogstad (1986). This method yields a unimodal or bimodal 
8(0;^ that exactly fits the measured buoy cross-spectra. Directional resolution is limited 
because the buoy cross-spectra define only four moments of The buoy spectral 

analysis is based on an ensemble average of 54 data segment with lengths of 200 s. The 
resulting E(^0) estimates with 0.005 Hz frequency resolution were smoothed to the same 
resolution as the array estimates. At low frequencies where they overlap, the array and 
buoy estimates are in good agreement. The final estimates of E(f,d) in 20 m depth were 
formed by matching the array estimates in the frequency range 0.04 - 0.2 Hz with the 
buoy estimates in the frequency range 0.2 - 0.5 Hz. 
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m. OBSERVATIONS 



Example observations of the shoaling evolution of EC/'jO) are shown in Figures 4 
- 1 1. The foiu cases described in detail below span a wide range of conditions, including 
normal and large obhque incidence angles and both non-brealdng and breaking conditions. 
In order to distinguish linear refraction effects from nonlinear harmonic generation and 
surf zone effects, the observed evolution of E(/‘,0) was compared to simple linear theory 
predictions. On a beach with straight and parallel depth contours, the transformation of 
EO^,0) is governed by 



0) 

5 gs 

(Longuet-Higgins, 1957; LeMehaute and Wang, 1982), where subscripts d and s denote 
wave properties at a deeper and shallower cross-shore location, C and Cg represent the 
wave phase and group speed, and the angles 0^ and 0j are given by Snell’s Law 



C 

sin(0p =— sin(0p 

d 



Eqs (1) and (2) were used to transform the 20 m depth estimates of E(/',0) to each of the 
shallower array locations. The resulting predictions are included in Figures 4-11. 
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A. 



AUGUST 17 



The August 17 conditions are characterized by a steady oblique-offshore wind with 
a speed of 8 m/s and direction of -140°. Significant wave heights were less than 0.5 m 
at all array locations and breaking occurred in the region shoreward of the shallowest 
array. Frequency-directional wave spectra observed at the six arrays are compared to the 
linear refraction predictions in Figure 4. Observed and predicted directional distributions 
of energy at selected frequencies are compared in Figure 5. In 20 m depth, the wave 
spectrum in the dominant 0.15 - 0.30 Hz range is broad in frequency and directionally 
narrow with waves propagating at large southerly incidence angles nearly parallel to the 
shore. These waves are strongly refracted over the inner shelf and beach, and the 
observed changes in wave direction are well described by the linear refraction predictions. 
High frequency spectral energy levels are reduced by about a factor 5 between the deepest 
and shallowest sites, consistent with predicted spectral levels (left panels). The observed 
refraction of the smaller and broader 0. 1 Hz swell peak from -50° at the deepest site to 
-20° at the shallowest site is also in agreement with the prediction (Figure 5, left panels). 
This peak is amplified during shoaling and dominates the wave spectrum at the shallower 
sites (Figure 4), in agreement with the prediction. Whereas observed energy levels and 
mean propagation directions are in good agreement with the predictions, the observed 
5(0;^ are broader than the predicted extremely narrow S(0/) (e.g., 0.17 Hz in Figure 5). 
Possible explanations for these discrepancies are the limited directional resolution of the 
array and alongshore variations in wave directions caused by small alongshore depth 
variations that might broaden the array estimates of 5(0;/). 
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B. 



OCTOBER 2 



The shoaling evolution of a more energetic wave field was observed on October 
2 and is detailed in Figures 6 and 7. Conditions are characterized by a steady onshore 
wind with a speed of 1 1 m/s and a direction of 55'. In 20 m depth, E(/j0) is relatively 
narrow in both frequency and direction with a single peak at approximately / = 0.17 Hz 
that is closely aligned with the wind direction (bottom panels of Figures 6 and 7). At 
higher frequencies, the S(Q'J) estimates (extracted from the directional wave buoy) are 
bimodal and approximately synunetric about the wind direction. The separation of the 
two peaks increases with frequency, similar to recent observations by Young et al. (1995) 
and Ewans (1998). Young et al. (1995) show that the bimodal stmcture of E(f,Q) at high 
frequencies can be explained with theoretical predictions of transfers of energy to higher 
frequencies through nonlinear wave-wave interactions. However, the 8 m depth array 
estimate of E(f,0) shows high frequency waves propagating in the local wind direction 
and suggest that the bimodal distributions estimated from the buoy data may be erroneous. 
Errors in the directional buoy estimates of 8(0;/) at high frequencies may be introduced 
through nonlinearities in the buoy response or the inherent lack of resolution of 
directional buoys that measure only four integral moments of 8(0;/) (e.g., Herbers and 
Guza, 1990). 

The observed significant wave height at the shallowest array is about 10% smaller 
than the predicted value, indicating that the surf zone is confined to a region shoreward 
of the shallowest array (upper left panel in Figure 6). The dominant sea peak is strongly 
refracted over the inner shelf and beach, and the observed changes in wave direction are 
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well described by the linear refraction predictions. Similar to the previous case, the 
predicted refracted peaks are directionally narrower than the observed peaks (left panels 
in Figure 7). At two times the spectral peak frequency (0.34 Hz) the E(f,0) observed at 
the shallower arrays are bimodal with one peak (0 « 50') aligned with the wind direction 
and the expected harmonic peak (0 = 20') aligned with the primary 0.17 Hz sea peak 
(Figure 6). This harmonic peak is strongly amplified during shoaling and dominates the 
high-frequency spectrum in 2.9 m depth (Figure 6). 

C. NOVEMBER 1 

The November 1 case (Figures 8 and 9) is characterized by a steady oblique 
offshore wind with a speed of 7 m/s and direction of -140', similar to the August 17 case. 
Observed significant wave heights ranging from 1.2 m at the deepest array to 1.1 m at 
the shallowest array are close to the linear theory predictions (left panels in Figure 8) 
indicating that breaking occurred in the region shoreward of the shallowest location, fri 
20 m depth, the E(f,0) estimate features a dominant 0.10 Hz peak with a large southerly 
incidence angle (-60'). The high frequency part of the spectrum is directionally narrower, 
with even larger, almost shore-parallel southerly angles (-80'), similar to the August 17 
case. The strong refraction over the inner shelf and beach of the dominant wave direction 
is again well described by the linear refraction predictions (Figure 9, left panels). As in 
the previous cases, predicted directional distributions are narrower than observed at the 
shallower arrays. At high frequencies, second (0.2 Hz) and third (0.3 Hz) harmonic peaks 
are observed at the shallower sites which are approximately directionally aligned with the 
dominant (0.1 Hz) waves (Figure 8). Similar to the October 2 case, these harmonic peaks 
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are amplified between 5.5 and 2.8 m depths and dominate the directional distributions at 
high frequencies at the shallowest array (Figure 9, center and right panels). 

D. OCTOBER 19 

The most energetic waves in the SandyDuck experiment were observed on October 
19 (Figures 10 and 11) during a severe Nor’easter with maximum sustained wind speeds 
of about 17 m/s. The surf zone extended several kilometers offshore with a significant 
wave height at the 20 m depth array of about 3.9 m, decreasing steadily to 1.56 m at the 
shallowest array. Propagation directions of the dominant 0. 1 Hz waves are nearly normal 
to the beach in 20 m depth and virtually imchanged throughout the surf zone (Figure 11, 
left panels) even though spectral energy levels are reduced by more than an order of 
magnitude at the shallowest array (Figure 10, left panels). The directional distributions 
are slightly broader than predicted with linear refraction theory at the shallower arrays, 
but the broadening of S(0^ is much less pronounced than that observed by Herbers et 
al. (1998) for comparably energetic waves breaking on a shallow bar. At twice the 
spectral peak frequency (0.19 Hz), the observed spectra show a transition from wind seas 
in 20 m depth (propagating in the local wind direction 0 « 80*) to second harmonic 
waves that are ahgned with the dominant wave direction (0 « 0°). In 8.2 m depth the sea 
and harmonic peaks are approximately equal in magnitude, and at the shallowest arrays 
the sea waves are completely submerged in the harmonic energy levels. A similar 
transition from wind-generated seas to harmonic waves is observed at 0.29 Hz (Figure 11, 
right panels). As in the August 17 case, the bimodal S(0i/) observed in 21.2 m depth 
may be a directional buoy artifact rather than a true feature of the wave field. In 8.2 m 
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depth, a refracting sea peak from 40* co-exists with a developing third harmonic peak at 
5° (ahgned with the dominant waves). This harmonic peak broadens during shoaling and 
dominates the spectrum at the shallower arrays. 
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IV. SUMMARY 



The evolution of the frequency-directional wave spectrum, E(/j0), across the inner 
continental shelf and beach was examined with extensive measurements collected at the 
Army Corps of Engineer’s Field Research Facility near Duck, NC during the recent 
SandyDuck experiment. The barrier island field site is characterized by a straight 
coastline with a gently sloping (1:250) sandy bottom and weak alongshore depth 
variations. Arrays of bottom pressure sensors were deployed on the inner shelf in 20 m 
depth (5 km from shore) and on the beach in depths ranging from 2 - 5 m (Figure 1). 
A directional wave buoy was added to the 20 m depth array to measure high frequency 
seas that are attenuated at the sea bed. The SandyDuck arrays were complemented by a 
permanent array of pressure sensors in 8 m depth maintained by the Army Corps of 
Engineers. High-quality wave data were collected nearly continuously at all instrument 
locations during a four month period (August - November 1997). 

The preliminary analysis presented here is focused on four case studies that span 
a range of conditions and illustrate the observed wave shoaling evolution in both non- 
breaking and breaking conditions. Smooth estimates of E(/,0) were estimated from the 
cross-spectra of six arrays using the variational technique described in Berbers and Guza 
(1990). In 20 m depth, the array estimates of E(f,0) at low frequencies were combined 
with buoy estimates at high frequencies obtained with the Maximum Entropy Method 
described in Lygre and Krogstad (1986). In all four case studies, the cross-shore 
evolution of the dominant wave propagation direction agrees well with simple linear 
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theory predictions for wave refraction on a beach with straight and parallel depth 
contours. In the more energetic cases, the observed E(f,0) show the development of 
harmonic peaks at the shallower sites with directions that are approximately aligned with 
the dominant waves (Figures 6, 8, 10), qualitatively consistent with theoretical wave- wave 
interaction rules and earlier observational studies (Freilich et al. 1990; Berbers et al. 
1998). 

In both non-breaking and breaking conditions, the observed E(/',0) are directionally 
broader than the linear theory predictions, in particular when the incident wave field is 
directionally narrow (e.g., August 17, Figures 4, 5). These discrepancies may reflect the 
limited resolution of the arrays and/or alongshore variations in wave propagation 
directions associated with weak alongshore depth variations. Whereas Berbers et al. 
(1998) observed a sudden, dramatic directional broadening of wave spectra when waves 
broke on a shallow submerged sandbar, the present observations on a nearly planar beach 
do not indicate a strong wave breaking effect on directional wave properties. In fact, in 
the most energetic case (October 19) when the surf zone extended across the entire 
instrumented transect, the directional width of E(/',0) remains nearly constant upon 
shoaling. These differences are possibly related to differences in nonlinear shoaling 
evolution and wave breaking characteristics (e.g., strongly localized breaking on a sand 
bar versus gradual breaking on a planar beach). Further work is needed to clarify the role 
of the sea bed profile in the shoaling evolution of directional wave properties. 
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APPENDIX 



Figure 1. Plan view of the experiment site with locations of the nearshore, 8 m, 
and 20 m arrays and associated bathymetry (depth contours are relative to mean sea 
level). 

Figure 2. Detailed plan view of the nearshore, 8 m, and 20 m arrays. Pressure 
sensors locations are represented by solid dots and the directional wave buoy location is 
indicated with a circle. 

Figure 3. Representative beach topography during SandyDuck (courtesy of the 
Army Corps of Engineer’s Field Research Facility). The bathymetry is characterized by 
a small shore-parallel bar located 275 m offshore and nearly straight and parallel depth 
contours. Bathymetry changes over the course of the experiment were weak with the 
exception of the dynamic beach face region shoreward of the nearshore arrays. 

Figure 4. Wave shoaling evolution from 19.8 m depth (bottom panels) to 2.5 m 
depth (top panels) of small amplitude waves arriving from large oblique southerly angles 
on August 17, 1997. Observed (blue) and predicted (red) frequency spectra and 
significant wave heights are shown in the left panels. The color contour panels show the 
observed and predicted frequency-directional spectra at each array (units cm^Hz"; 0=0° 
corresponds to normal incidence to the beach) with the water depths indicated in the 
upper right comer . 

Figure 5. Observed (soUd) and predicted (dotted) directional distribution of wave 
energy at selected frequencies on August 17, 1997. Predictions at frequencies higher than 
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0.2 Hz (only measured with the directional wave buoy in 20 m depth) are omitted. 

Figure 6. Wave shoaling evolution from 20.3 m depth (bottom panels) to 2.9 m 
depth (top panels) of moderately energetic waves arriving from large oblique northerly 
angles on October 2, 1997 (same format as Figure 4). The frequency-directional values 

o 

are multiplied by / to enhance the directional properties of high frequency waves. 

Figure 7. Observed (sohd) and predicted (dotted) directional distribution of wave 
energy at selected frequencies on October 2, 1997 (same format as Figure 5). The local 
wind direction is indicated with a dashed vertical line. 

Figure 8. Wave shoaling evolution from 20.8 m depth (bottom panels) to 2.8 m 
depth (top panels) of moderately energetic waves arriving from obhque southerly angles 
on November 1, 1997 (same format as Figure 4). The frequency-directional values are 
multiplied by / to enhance the directional properties of high frequency waves. 

Figure 9. Observed (soUd) and predicted (dotted) directional distribution of wave 
energy at selected frequencies on November 1, 1997 (same format as Figure 5). 

Figure 10. Wave shoaling evolution from 21.2 m depth (bottom panels) to 3.3 
m depth (top panels) of large amplitude waves arriving from directions close to normal 
incidence on October 19, 1997 (same format as Figure 4). 

Figure 11. Observed (solid) and predicted (dotted) directional distribution of 
wave energy at selected frequencies on October 19, 1997 (same format as Figure 5). The 
local wind direction is indicated with a dashed vertical line. 
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